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Motivated by recent experiments with Bechgaard salts, we investigate the competition between 
antiferromagnetism and triplet superconductivity in quasi one-dimensional electron systems. We 
unify the two orders in an SO(4) symmetric framework, and demonstrate the existence of such 
symmetry in one-dimensional Luttinger liquids. SO (4) symmetry, which strongly constrains the 
phase diagram, can explain coexistence regions between antiferromagnetic, superconducting, and 
normal phases, as observed in (TMTSF)2PF6. We predict a sharp neutron scattering resonance in 
superconducting samples. 



A common feature of many strongly correlated elec- 
tron systems is proximity of a superconducting state to 
some kind of magnetically ordered insulating state. Ex- 
amples include organic materials heavy fermion 
superconductors 0,01, and high Tc cupratesly. Several 
theoretical analyses suggest that a strong repulsion be- 
tween the two orders plays an important role in determin- 
ing the phase diagram and low energy properties of these 
materials 0. The idea of competing orders was devel- 
oped into the SO (5) theory of high Tc superconductivity 
of S.C. ZhangQ. S0(5) symmetry has also been apphed 
to study competition of ferromagnetism and p-wave su- 
perconductivity in Sr2Ru04@|, and antiferromagnetism 
and d-wave superconductivity in k-BEDT-TTF salts 0. 

In this paper we consider the interplay of antiferro- 
magnetism (AF) and triplet superconductivity (TSC) 
in quasi one-dimensional (QID) electron systems. Our 
study is motivated by QID Bechgaard salts (TMTSF)2X. 
The most well studied material from this family, 
(TMTSF)2PF6, is an antiferromagnetic insulator at am- 
bient pressure and a superconductor at high pressures [lol 
IllL I12L [13 . The symmetry of the superconducting order 
parameter in (TMTSF)2PF6 is not yet fully established, 
but there is strong evidence that electron pairing is spin 
triplet: the superconducting Tc is strongly suppressed by 
disorder|l4l|: critical magnetic field Hc2 in the interchain 
direction exceeds the paramagnetic limit [3] ; the electron 
spin susceptibility, obtained from Knight shift measure- 
ments, does not decrease below Tc[3- I^i another mate- 
rial from this family, (TMTSF)2C104, superconductivity 
is stable at ambient pressure and also shows signatures 
of triplet pairing [13 lla, lid l2d|- Insulator to supercon- 
ductor transition as a function of pressure has also been 
found in(TMTSF)2AsF6|2ll. 

The phase diagram of interacting electrons in one di- 
mension was obtained in Ref. 113 using bosonization and 
renormalization group (RG) analyses. At incommensu- 
rate filling, this system has a phase boundary between 
spin density wave (SDW) and TSC phases when Kp = 1 
and gi > {Kp is the Luttinger parameter in the charge 
sector, gi the backward scattering amplitude). The start- 
ing point of our discussion is the observation that, in 
the absence of umklapp, ID Luttinger liquids have an 



"isospin" S0(4)iso symmetry |23| at the boundary be- 
tween SDW and TSC phases. To define this symmetry, 
we introduce the charge of left and right moving elec- 
trons, Q±, and two new operators, 0±, (r = ±) 



Qr 



-kl- 



(1) 



Here Oj. creates right/left moving electrons of momen- 
tum ±kf + k and spin s. Combining these according to 

= (ot + e,)/2, = (et - e,)/2i, and j; = q„ 

we see that the generators satisfy two independent chi- 
ral SO (3) algebras, \ja,Jb] ^ i5'''''' e''^" J^. The product 
of left and right algebras yields the total isospin group 
S0(4)iso wSO(3)flxSO(3)L. The Luttinger Hamiltonian 
at incommensurate filling, 7i, generically commutes with 
the charge operators Q±. Using a bosonized form of TL, 
it can be shown that at Kp — 1, the 0± operators also 
commute with Thus, at Kp = 1, S0(4)iso 

forms an exact symmetry of 7i. In addition, for spin- 
symmetric interactions, which describe Bechgaard salts 
to a very good approximation |3 the system has 



S0(3)spin symmetry, [Sa 
the total spin 



'/3J 



^^'^Sj, generated by 
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Hence, for Kp — 1, i.e. the line separating SDW and TSC 
phases, the system has full SO(3)spinXSO(4)iso symme- 
try. We emphasize that for Luttinger liquids at incom- 
mensurate filling, this symmetry always appears at the 
SDW/TSC phase boundary and does not require fine- 
tuning of the parameters. 

This symmetry can be used to unify SDW and TSC 
order parameters. SDW order away from half filling is 
described by a complex vector order parameter. 



(3) 
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QID band structure restricts the orbital component of 
triplet superconducting order to be ^'(p) oc p^;, with x 



2 



the coordinate along the chains. Thus, the TSC order 
parameter is also a complex vector, 

- ^E«l.fe.(^"^2)..'aL__,,, (4) 

kss' 

The four vector order parameters Re$, Im<i>, Re^, and 
Im\l/ can be combined into a 4x3 matrix, 

/(Re*), (Im*), (Re$), (Im$), \ 
Paa = (Re#)y (Im#)j, (Re$)y (Im*)^ (5) 
V(Re^), (Im^), (Re$), (Im$), / 

Each column (row) of P transforms independently as a 
vector under the action of S0(3)spin (S0(4)iso)- 

The Qr operators are reminiscent of Yang's rj opera- 
tor, which generates an S0(4) symmetry for the Hubbard 
model |2^. Unlike rj, whose center of mass momentum is 
always the commensurate wave vector it, the @r have 
their momenta at the wave vectors ±2fc/. This is crucial 
for defining the symmetry at arbitrary electron density; 
in contrast, Yang's SO (4) applies only at half-fihing. 

Real materials are only QID and coupling between 
chains gives rise to finite temperature phase transitions. 
However, as long as 3D coupling is weaker than the in- 
trachain tunnelling and interactions, the nature of the 
ordered state is determined by the most divergent suscep- 
tibility within individual chains. Hence, for QID materi- 
als near the SDW/TSC boundary we expect to find Kp 
close to one, and to find approximate S0(3)spin xS0(4)iso 
symmetry. Then, the phase diagram in three spatial di- 
mensions is obtained from a Ginzburg-Landau (GL) free 
energy whose form is strongly constrained by symmetry, 

F = ^{VP-a^f + fPi^ + 5r (P2^ + Pi - - Pi) 

+ U,PlPlp+U2P-ac.P-apPlo.Plp (6) 

This is the most general expression with S0(3)spinX 
S0(4)iso symmetric quartic coefficients. We follow the 
common assumption that changing the external control 
parameters of the system, e.g. temperature and pres- 
sure, only affects the quadratic coefficients. These are 
thus allowed to break the symmetry and tune the phase 
transition. For 5r ^ 0, the symmetry is broken down to 
SO(3)spi„xSO(2)cXSO(2)t, where S0(2)c is generated by 
the total charge Q = Q+ -I- Q-j and S0(2)t comes from 
the lattice translational svmmetrv|2^. 

At half-filling, umklapp scattering must be taken into 
consideration. This is the case for Bechgaard salts, where 
structural dimerization splits the conduction band into a 
full lower band, and a half-filled upper band. Umklapp 
turns two right movers into left movers, and vice versa, 

'^3 ||H"+,fc+gs«+,p-9ta-,p*a-.fcs + /I.C. (7) 

The term ((Jj) pins the phase of the SDW, reducing it 
to the real Neel vector N = Re^ (for > 0). In the 



presence of umklapp there is still a direct AF/TSC tran- 
sition at Kp = ljl|. On the other hand, Q+ and Q- 
are no longer conserved separately, and S0(4)iso sym- 
metry of the free energy is broken. However, to linear 
order in g^, S0(4)iso is not broken all the way down to 
S0(2)c. The contribution of Ti,3 to the free energy is 

AF = ^ ^(Re$)^ — (Im$)^^ , which preserves a resid- 
ual symmetry S0(3)iso, given by the diagonal subgroup 
of S0(4)iso «SO(3)j?xSO(3)l|23. The generators of 
S0(3)iso are = J+ + J". Together with S0(3),pi„ 
invariance, the total symmetry defined by these opera- 
tors is S0(4)«S0(3)spi„xS0(3)iso. 

The justification for considering small g^ is as follows. 
The bare value of 53 is proportional to dimerization, 
which is only of order of 1% in (TMTSF)2PF6^]. Fur- 
thermore, the GL free energy depends on the renormal- 
ized value 53, off at the crossover scale between ID and 3D 
physics, which is even smaller than the bare value of (73. 
Umklapp is irrelevant inside the TSC phase, as well as 
on the AF/TSC phase boundary. Even in the AF phase, 
where umklapp is a relevant perturbation, g^ fiows near 
zero before diverging, and this divergence may be cut 
off by the onset of 3D coupling. Therefore, everywhere 
close to the AF/TSC phase boundary, we can assume 
that (73, cff is small. 

SO (4) symmetry unifies AF and TSC orders, which are 
now combined into a 3 x 3 tensor order parameter, 

/ (Re*), (Im*), \ 
Qaa = {Re% (lm*)y Ny (8) 
V (Re*), (Im*), N, J 

The columns (rows) of Q transform as a vector under 
the spin (isospin) SO (3) algebra, [Sa,Qbi3] = ie^'^'^Qb-y 
i[Ia,Qbf3] — ie°''"'Qci3)- In analogy with ©, the GL free 
energy near the AF/TSC phase boundary is 

+ ■UlQlaQbl3 + QapQbaQbp (9) 

When 5r — Q the model has full SO (4) symme- 
try. Away from this line it only has spin and charge 
SO(3)spinXSO(2)c symmetry. A derivation of the GL 
free energy for weakly interacting QID electrons yields 
the model in © with ui = 2lC,{^)/lQiT'^VfT'^ and U2 = 
-7C(3)/87r2v^r2:24]. 

The properties of model depend strongly on the 
sign of U2, which determines whether the triplet su- 
perconductor is unitary {^ip oc Si/") or non-unitary 
(3?^ X 3?/' / 0). We expect the unitary case, < 0, 
to be of experimental relevance to (TMTSF)2PF6 , and 
in the remainder of this paper we primarily concentrate 
on it. The mean field diagram is then composed of an 
AF phase separated from a TSC phase by a first order 
phase transition, and a disordered (Normal) phase sep- 
arated from the two other phases by second order lines 
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FIG. 1: Mean field phase diagram of eq. @ in the unitary case 
■112 < 0. There is a first order transition (thick fine) between 
AF and TSC phases. Inset: Corresponding diagram for non- 
unitary case U2 > 0. M denotes a mixed AF/TSC phase. 




FIG. 2: Renormahzation group flow of the SO(4) symmetric 
theory eq. Q in d = 4 — e. There are no stable fixed points. 
Instead, there are two types of runaway flow, corresponding 
to unitary (u-z < 0) and non-unitary (w2 > 0) TSC. 

(see Fig. QJ. For completeness, we include in the inset 
the mean field phase diagram for the non-unitary case. 

To understand the role of thermal fluctuations in 
model Q and in slightly perturbed models where the 
quartic coefficients do not lie exactly on the SO (4) sym- 
metric manifold, we use 4 — e RG analysis. We find 
that the RG equations have only two fixed points: a 
trivial Gaussian fixed point, f ^ 5r = ui = 0, and 
an S0(9) Heisenberg point, f 7^ 0, (5r = 0, ui 7^ 0, 
U2 = 0. In Fig. El we show RG flows in the SO (4) 
symmetric plane, where we find runaway flows whenever 
"5,2 7^ 0. The analysis can be generalized to order param- 
eters and ^ that are A^-component vectors, in which 
case the SO(4)«SO(3)spinXSO(3)iso symmetry becomes 
S0(A^)spi„xS0(3)iso. We find that even in the large N 
Hmit, ah flows with U2 < are runaway flows, indicating 
the absence of fixed points with unitary TSC. 

The absence of a fixed point in the RG flow often im- 
phes that fluctuations induce a first order phase transi- 
tion, thus precluding a multicritical point in the phase 
diagram. In order to inspect this possibility, we study 
model (jnj directly in 3D in the large N hmit Hi]. The 
idea of the large A^ expansion is to sum self-consistently 
all bubble diagrams. We find (inset of Fig. |3Jl a first 
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FIG. 3: Schematic temperature-pressure phase diagram of 
(TMTSF)2PF6 El El. AF/N and AF/TSC correspond to 
coexistence regimes of the appropriate phases. Inset: Phase 
diagram for competing AF and unitary TSC states for model 
@ in the large N limit. 5r tunes the system across the two 
phases. Thick lines represent first order transitions. 



order transition between AF and TSC phases along the 
SO (4) symmetric line 5r = 0, in agreement with mean 
field theory. A new feature of the large A^ limit is the 
first order transition between the Normal and AF phases 
in the vicinity of the critical point. If we assume that the 
experimentally controlled pressure changes an extensive 
variable conjugate to 5r, such as the volume of the sys- 
tem, the first order transition broadens into a coexistence 
region of TSC and AF. This is consistent with the exper- 
imental phase diagram of (TMTSF)2PF6 (Fig. ^. An 
unusual feature of the theoretical phase diagram (inset. 
Fig. O is the first order transition between Normal and 
TSC phases. This is similar to the fluctuation-driven first 
order transition between Normal and superfluid phases 
proposed for ^He by Bailin et al^^. In bulk ■^He, the 
coherence length is very long and the transition is mean- 
field like. The fluctuation region is so small that tiny 
discontinuities at the transition caused by fluctuations 
would be impossible to observe. In QID systems, such 
as Bechgaard salts, the fluctuation region is expected to 
be large[30| and the discontinuous nature of the Normal 
to TSC transition may be experimentally accessible. This 
transition has been investigated through specific heat 
measurements in (TMTSF)2C104 3l|. The results were 
interpreted as a mean-field BCS transition, although the 
amplitude of the specific heat jump was unusually large. 
Since in these materials is very sensitive to impurities, 
the extra amplitude in the specific heat jump might be 
attributed to broadening due to disorder of a (5-function 
peak in specific heat. Experimental observation of the 
coexistence region of Normal and TSC phases, or of hys- 
teresis effects in resistivity measurements may also be 
very difficult. Nearly equal strains in the TSC and the 
Normal phases and the nearly temperature independent 
superconducting transition temperature can make such 
a coexistence region very small. On the other hand, it 
is possible that spin anisotropy present in real materi- 
als, but not included in our theoretical analysis, stops 



4 



the runaway RG flows leading to the first order tran- 
sition. Another possibihty is that one loop 4 — e RG 
calculations and the large N expansion do not capture 
the correct behavior of model © for e = 1 and = 3. 
For the Normal to unitary TSC transition, De Prato et 
al. argued that a stable fixed point describing the second 
order phase transition appears in the six-loop expansion 
of the GL free energy[23- We hope that future experi- 
ments will investigate the nature of the Normal to TSC 
transition in Bechgaard salts in more detail. 

The most dramatic consequence of enhanced symmetry 
at the phase transition is the prediction of new low en- 
ergy collective excitations. In the vicinity of the AF/TSC 
phase transition, SO (4) symmetry leads to a new collec- 
tive mode, corresponding to the 6 operator that rotates 
AF and TSC orders into each other [2J|. As pressure is 
varied toward the phase transition, breaking of the SO (4) 
symmetry is reduced, leading to a decrease in the en- 
ergy of the 8 excitation. Mode softening is not expected 
generically at a first order phase transition and identifies 
the 0-resonance as a generator of the SO (4) quantum 
symmetry. Weak symmetry breaking due to interchain 
coupling and higher order umklapp terms may lead to 
a small gap and to finite broadening of O, even at the 
AF/TSC phase boundary. Deep in the normal phase the 
Q excitation cannot be probed by conventional methods, 
such as electromagnetic waves or neutron scattering, as 
these only couple to particle-hole channels (e.g. spin or 
density) and O is a collective mode in the particle-particle 
channel. The situation changes when the system be- 
comes superconducting. In the presence of a condensate 
of Cooper pairs, charge is not a good quantum number 
and particle-particle and particle-hole channels mix. The 
O excitation should thus appear as a resonance in inelas- 
tic neutron scattering experiments |24| . and its intensity 
should be proportional to the square of the superconduct- 
ing amplitude l^'p. For QID Bechgaard salts, we expect 
strong pairing fluctuations even above Tc. Hence, pre- 
cursors of the O resonance may be visible in the normal 
state, with strong enhancement of the resonant scattering 
intensity appearing when long range TSC order develops. 

It is useful to put the SO (4) model of AF/TSC com- 
petition in Bechgaard salts in the general perspective 
of electron systems with competing orders. In the 
case of the SO (5) theory of AF and d-wave SC in 2d 
systems 0,0 1 it is difficult to construct realistic micro- 
scopic models with such symmetry (see e.g. Refs.|3IH)- 
By contrast, S0(4) symmetry in Bechgaard salts arises 
naturally from a conventional Luttinger description. We 
also point out that tuning across the AF/TSC phase 
boundary in (TMTSF)2PF6 can be done in the same 
sample by varying pressure, whereas tuning the AF/SC 
transition in the cuprates requires using different sam- 
ples. Thus, we consider (TMTSF)2PF6 a good candi- 
date for experimental observation of emergence of higher 
symmetry in a strongly correlated electron system. 



In the discussion above, we assumed Luttinger liquid 
behavior in individual chains to motivate the approx- 
imate SO (4) symmetry at the AF/TSC phase bound- 
ary. It has been suggested that for the superconducting 
phase of Bechgaard salts, interchain hopping is strong 
enough to turn the system into a strongly anisotropic 
Fermi liquid [s^l- The decreased nesting condition in 
this case strongly affects antiferromagnetism. For the 
classical symmetry of the GL free energy, this effect can 
be absorbed into the normalization of the field N, so 
that the GL parameters only display a small deviation 
from SO (4) symmetry at the mean-field level. Thus, 
we do not expect a qualitative change in the phase di- 
agram presented in Fig. |3| (see Ref. for a detailed 
discussion). To verify the approximate quantum SO (4) 
symmetry for the strongly anisotropic Fermi liquid, one 
can study the spectrum of collective excitations using 
an RPA-t ype analysis and verify the existence of the O 
excitation|38j. These results will be presented elsewhere. 

In summary, we introduced an SO (4) framework for 
the competition between AF and TSC in QID electron 
systems. The microscopic origin of the SO (4) symme- 
try at the transition between AF and TSC orders was 
identified in the Luttinger liquid model. Our results 
have direct implications for QID organic superconductors 
from the (TMTSF)2X family. For example, first order 
transitions between AF and TSC phases, and between 
AF and Normal phases, explain the AF/TSC and the 
AF/Normal coexistence regions found in the phase dia- 
gram of (TMTSF)2PF6 ,12]. We also argue that the Nor- 
mal/TSC transition in these materials could be weakly 
first order. We predict a sharp resonance in neutron scat- 
tering, whose characteristics identify it unambiguously as 
a generator of SO (4) symmetry. 
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